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We propose a model for the antihyperon polarization in high-energy proton-nucleus inclusive re- 
actions, based on the final-state interactions between the antihyperons and other produced particles 
(predominantly pions). To formulate this idea, we use the previously obtained low-energy pion- 
(anti-)hyperon interaction using effective chiral Lagrangians, and a hydrodynamic parametrization 
of the background matter, which expands and decouples at a certain freezeout temperature 
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^ ; I. INTRODUCTION 

The unexpected discovery, by Bunce et al. in 1976 of large A polarization in the inclusive production, induced by 
proton-Beryllium collisions at 200 GeV, opened an interesting field of investigation. Earlier experimental data, as those 
on the TTp, pp and Kp collisions Q, indicated that in general the polarization decreased with the energy and, moreover, 
theoretically it was hard to believe that high-energy inclusive reactions could produce significant polarization, because 
, it was expected that a large cancellation occurred when many uncorrelated channels are superposed. 
■ In the following years many proton-nucleus collision experiments have been performed and these experiments 
I ' confirmed the A polarization, in the pt range of 0-2 GeV 0, Q . They also showed that in the same domain the A 
polarization is very small, consistent with 0. This fact created the false idea that only the hyperons could be polarized 
[ — ' and the antihyperons not. In 1981, Wilkinson and his collaborators measured for the first time the polarization of 
, a charged hyperon, S"*" |^, and found to be positive. In 1983 0, the S~ polarization was measured, with results 
'n ■ consistent with the S+ ones, and in 1986, the S~ polarization was observed, and it was found to be negative 
, More recently, polarizations of hyperons produced with a E~ beam also have been measured 

Soon it was possible to obtain polarizations of other hyperons: in 1990 the S was measured negative, and 

in 1993 Morelos and his collaborators measured the S polarization that is positive. These results showed that, 
Oh refuting the earlier belief, the antihyperon polarizations are not always zero. 

^ Many models have been proposed trying to explain these experimental res ults, such as the Lund model [l^, the 

quark recombination model 13], and the OPER (reggeized one pion exchange) [TJIj that explained the A polarization 
but were not able to calculate the A polarization. These models are based on the leading-particle effect, that means 
that the observed hyperon is considered as a direct product of the incoming proton after a recombination. Clearly, 



this mechanism can be effective only for hyperons but not for anti-hyperons. 
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S ■ Later, models based on constituent quarks 'l5l and others ^,1121 were proposed, but a satisfactory explanation 
for the antihyperon polarization was not achieved. As we can see, the models are able to explain the A polarization 
(and even A), but in fact, all of the models predict the production of unpolarized antihyperons and this fact shows 
that these models must be supplemented in order to describe also the antihyperon data. 

It is well known that final-state interactions are fundamental in the study of the nonleptonic hyperon decays, 
where the final-state strong phases must be considered 0|, |0, 0|. In high-energy inclusive processes, a large 
number of particles is produced, so we expect that something similar occur in this kind of reactions. In 1993, Y. 
Hama and T. Kodama used the hydrodynamical model to study the effect of final-state interactions of the 
produced (anti) hyperons with the background matter. It was supposed that the (anti) hyperons are produced inside 
a medium composed of hot hadronic matter, that is formed during the collision, and the polarization is the result of 
the interaction of the hyperon with the surrounding matter. In [2l|, this interaction was represented by an optical 
potential and they showed that this mechanism could produce polarized hyperons and antihyperons. However, some 
important questions remained opened, namely, the origin of the potentials and why different potentials are needed 
for the different (anti) hyperons. 

The purpose of this paper is to improve the model described in the preceding paragraph, by explicitly including some 
detailed microscopic interaction between the antihyperon and the surrounding thermalized medium J,]. Considering 
that in this kind of collisions pions are predominant among the produced particles, as a first approximation, we will 
include only the pion-antihyperon interactions in the study. The observed polarization is the average effect of these 
interactions. An important feature is that as the antihyperons and pions are in thermal equilibrium inside this fluid. 
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the energy involved in the microscopic interactions is generally low. Also, since the hyperons are produced mostly as 
the leading particles, we shall confine ourselves only to antihyperon polarization. 

A crucial factor to the polarization mechanism is the shape of the fluid rapidity distribution, dN/da, at the moment 
of emission. If dN/da were constant, no polarization would be expected because, for antihyperons with rapidity y, 
the contributions coming from the a > y region would cancel out the ones from the a < y domain, as they have 
opposite signs. But in the relevant a domain for the experimentally produced antihyperons, dN/da is decreasing, 
and then, the average polarization is dominated by the contributions from a < y. Another remark is that, since we 
are interested in antihyperons which are escaping, the microscopic collisions is essentially in the forward direction. 

The content of this work is the following. In the next section, the previously developed low energy pion-hyperon 
interactions [l^jlH,!^, based on chiral Lagrangians, are shortly reported. In Sec. Ill these results are used to calculate 
the polarizations in high-energy inclusive processes. In Sec. IV the results are shown and, finally, conclusions are 
drawn in Sec. V. 



II. LOW-ENERGY PION-HYPERON INTERACTIONS 



In this section, the basic question to be answered is how to describe the low energy t:Y interactions in a reliable 
way. The path to be followed is to study the ttY interactions and then to use the CPT invariance (neglecting the 
possible effects of CP violation) to determine the amplitudes in the ttY scattering. We will study the ttA, ttS and ttS 
interactions, and then, the quantities of interest that are cross sections and polarizations, will be calculated. The Y 
polarization, for example, is just the Y one, with the opposite sign. The cross sections for charge conjugate reactions 
are the same. 

Unfortunately, one must remark that this kind of interaction is not very well studied experimentally, the first 
measurement of phase shifts that has been made, was the difference Ss — Sp in the it A interaction for ^/s = m^, by 
the Fermilab E756 experiment and improved last year by the HyperCP experiment but still with reasonable 
experimental uncertainties. 

So, the option is to describe these interactions with a model. Some models [l^-ls^ have been proposed to study 
the ttA phase shifts, but these values are not accurately determined yet. On the other hand, the low energy ttN 
interactions can be described with a high degree of accuracy, with models based on chiral symmetry p^ . |28j |. So, 
in this work we will use a model ^J-_24J which is based on effective chiral Lagrangians, that considers baryons, spin 
3/2 resonances, rho mesons and the a term. An important feature of the model is the inclusion of resonances in the 
intermediate states. In the tt+P scattering (3/2 isospin channel) the A+"'"(1232) dominates the total cross section 
at low energies and is also very important in the other isospin channels. In the ttY case, we expect that the same 
behavior occurs. The Lagrangians to be considered are 
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where B, R^, 0, p are the baryon, the spin 3/2 resonance, the pion and the rho fields with masses tub, w^, /i, and 
rUp , respectively, /is is the baryon magnetic moment, are the isospin matrices and Z is a parameter representing 
the possibility of the off-shell-resonance having spin 1/2. 

The first case to be considered is the The ttA interaction. Since A has isospin 0, the scattering amplitude T^a has 
the general form 

T:^l^uip,)^^A+^-^^B'j5,MP) , (5) 

where and are the initial and final 4-momenta of A, and k'^ are those of the pion. Indices a and b indicate 
the initial and final isospin states of the pion. Fig. 1 shows the relevant diagrams, where we have omitted the crossed 
diagrams, although included in the calculations. To calculate the diagrams from Fig. 1, we must use the Lagrangians 
and |2Jl, that in this case are 
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/:a.s = ^ {S7m75tA} .9^0 + h.c. (6) 
2mA 



fA} .d''(l) + h.c. , (7) 
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FIG. 1: Diagrams for vrA Interaction 
Calculating the contribution from Fig. la (intermediate S), we have 
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The interaction with the intermediate resonance E*, shown in Fig. lb, is 



A^. = 4^|^i-!^l±^^(2m|.+mAms.-n.i + 2M^ 

+ [(mA+ms.)^+(2ms. + mA)Z^] /c.fc'l , 
Ss* = ^^f^ ^ , 9 



with 



A ^ - {2m|. - 2mi - 2mA™|. 
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2TO|ms«+Ai^(2mA-ms«)} + ^('7ia+"Ts*)* 
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B = ^„ [(m|, - m\)'^ ~ 2mAms-(TOs* + ^a)^ 
2m^ 
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(10) 



where u and are defined in the Appendix A. Note that eq.l^DJ shows small differences when compared with the 
one shown in |22| , where there was a mistake. The correct expression is presented here. 

In addition to these contributions,the a term (Fig. Ic) was also taken into account. In [s^, [s^, in the study of 
ttN interactions, it was included just as a parametrization 

Aa- = a + bt ^ 

B, = , (11) 

where a — 1.05 /i^^ and b — — 0.80/i^'^ are constants. Recent works |23|, js^ consider that the a term may be 
understood in terms of the exchange of two pions, in loop diagrams, as it is shown in Fig 2. 
So, these ideas lead to a cr contribution of the form 

Aa = crs + CTS* 

Ba = , (12) 
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FIG. 2: The scalar form factor receives contributions from tree interactions (white blob) and triangle diagrams with spin 1/2 
and 3/2 intermediate states. 

where a^, and cts* have been calculated in [23|. 

One must remark that considering the diagrams of Fig. 1 , the resulting amplitude will be real and consequently 
the S matrix will not be unitary. So, the unitarization of the amplitudes must be done. Observing that the scattering 
matrix may be expressed as 



, , Tba {a.k){a.ki) 

Mba = F = /l H TT- j2 



kkl 



(13) 



the partial wave decomposition can be made with 



ai± = 7; [Pi{x)fi{x) + Pi±i{x)f2{x)] 



(14) 



and the unitarization of these amplitudes, can be made if we reinterpret them as elements of K matrix and write 



ai± 



1 ~ ik an 



(15) 



The phase-shifts are then computed as 



\k au 



(16) 



and details on this calculation are shown in j22|, |2J|. 

The masses and magnetic moments that we used may be found in |36j |. The AttS coupling constant is g\-,^^=12.92 
(see Appendix A) . The resonance coupling constants were calculated with the comparison of the calculated amplitude 
and the Breit-Wigner one, with the data from 36] as it was done in 22] . The same procedure was adopted to calculate 
the other Yt:Y* couplings. The obtained value is qatt'S' is 9.38 GeV~^. With these results, we can calculate at, da/dfl 
and P, that will be needed in the next section. 
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FIG. 3: Total ttA cross section. 
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FIG. 4: Polarization in the ttA interaction, x —cos9. 



Observing Fig. O one can see a peak in the S*(1385) mass, what shows that the resonance dominates the total cross 
section at low energies. 

In the case of ttS interaction, the particles have isospin 1, so the composed system can have isospin 2,1 or 0. For 
this reason, the scattering amplitude in this case has the general form 

Taj^f3S = {TTyI]s\T\Tra^i3) 

= um{[A + il±lD.A']S^^d^s + [B + ^^^^^B']d^^df3s + [C+ ^^^^C"]S^56py}u{p) , (17) 
where a, /?, 7 and 6 are isospin indices of the particles. 
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FIG. 5: Diagrams to ttE Interaction 



The diagrams that we consider for the ttE interactions are shown in Fig. The interactions with the intermediate 
resonance I]*(1385) could also be included, but its decay branching ratio to the ttE channel is only 11%, what may 
be considered as a correction, and then we will neglect it. 

The amplitudes of eq. (|17() for the diagrams of Fig. [Slhave been calculated in 22], and the a contribution is given 

by 



A, = =a = = , (18) 

where a a, ctt. and a a* are taken from [23l |. 

The calculated total cross sections with E+ in the final states are shown in Fig. |H1 The resonance, that in this 
case is A*(1405), is still important, but now, the peak is not so high than it is in the other cases (less than 30 mb). 
It appears in the reactions where the isospin state is important. The polarizationsfor these reactions are shown in 
Fig. 13^ 

Now, let us turn our attention to the tt^ interaction. This case is very similar to the ttN scattering, because 
has isospin 1/2 (as the nuclcon) and the main difference is that the resonance of interest S*(1533) has isospin 7=1/2 
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(instead of 7=3/2 as A(1232)). Then, the scattering ampUtude has the general form 

n% = u{pi){[A++'^-t^B+]5,, 

+ [A- + ^-MlB-]ze,,,T^}u{p) . (19) 
The calculated and for the diagrams of Fig. |Slmay be found in p^.p^. 

In Fig. |51 we can see that one more time the resonance is very important, now the S(1533) resonance contribution 
dominates three of the reactions. 

So, now we have the ttA, ttS and ttS scattering amplitudes, that are the basic elements to calculate the polarization 
in high-energy processes. 



III. POLARIZATION IN INCLUSIVE PROCESSES 



In this section, it will be shown how to calculate the Y polarization in high energy processes. As it was said, we 
consider that during the collision, a system composed of expanding hot hadronic matter is formed. This system may 
be considered to be a fluid, composed of elements, that expands according to a Gaussian law |37j.|3^ 
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where a and at are the longitudinal and transverse rapidities of these elements of fluid, (f> is the azimutal angle and 
/3, and are parameters that determine the shape of this distribution. The values of these parameters depend on 
the participants of the collision (proton, beryllium, ...) and on their energies. 

The pions and hyperons are produced thermally inside these fluid elements, with initial momenta (relative to the 
fluid element) tTq (pions) and Aq (hyperons) and energies E', E'q, obeying the statistical distributions 
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The hyperon emerges with flnal momentum A', energy E' and polarization P' 
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FIG. 7: Polarizations in the ttE Interaction. 
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FIG. 8: Diagrams to tt^ Interaction 
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FIG. 9: TTi cross sections for channels with ^ in the final state. 



The produced hyperons interact, and as the relative energy between the tt and the Y is small, the amplitudes 
obtained in the last section may be used to calculate the polarization P' and da/dt. Then the average polarization 
may be calculated by the expression 



(23) 



that considers the elastic and with charge exchange Ri reactions, in the interactions of the hyperons with tt"*", tt" and 

The factor Q that appears in eq. (|23|l contains the statistical weighs of the production of the particles and the ones 
relative to the expansion of the fluid, and can be written as 



[d^N/da dat) 



(exp - l) (cxp (§) + 1 



^0 7^0 



^5[E', + E'-E'-dml + {7:', + K',-K'f] , 



and the eight dimension integration element is 



dr = da dat dAgdTTg 



(24) 



(25) 



In order to calculate the average polarization (|23|1 . wc must know the parameters /3, (3t and ag of the rapidity 
distribution (|20l) of the fluid elements and if this expression is a good choice. It is possible to answer these questions 
if we use the experimental data of the rapidity distributions of the pions produced in proton-proton collisions [s^ 
(that is a collision similar to the p— Be). In the hydrodynamical framework, the longitudinal rapidity distribution of 
the produced pions is given by 



d<T f da , da 
dy J da dy' 



(26) 
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FIG. 10: Polarizations in the ttH Interaction. 



that is a convolution of the rest-frame rapidity distribution da/dy', and the fluid rapidity distribution, da/da. 

The pion momentum distribution with respect to the fluid element (considering that they are produced thermally) 

is 

dp'^ E'^dy'df, ^£t_^ ' 

that after some manipulations gives 

- Ce-'''' (28) 

with p' ^ 0.98 for the temperature T ^ m^. 

The distribution da /da may be obtained from eq. H2l)|) . and will be a sum of Gaussian functions 

^ =^'[g-/3(o-ao)^ +g-/3(oW] ^ (29) 

da 
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FIG. 11: Parametrization of dN/dy for pp collisions compared with the data from [3fl| |. 

inserting this expression in H26|l . integrating and fitting the parameters A' , (3 and ao, we can compare the resuhs with 
the experimental data from ISR [s^ ■ The results are shown in Fig. 1111 

The best fits for da/dy are obtained with the parameters A' ~ 0.45, ao = 1-8 and (3 ~ 2.0 for 491 GeV 
GeV) and A' = 0.58, ao = 1.9 and f3 = 1.5 for 1030 GeV (y's7V7V=44 GeV) incident protons. Our real interest is for 
330, 400 and 800 GeV incident protons, that are the values for which there are antihyperon polarization data. The 
respective parameters may be obtained by interpolation. Observing Fig. ^]one can see that the /3 parameters may 
vary from 1.0 to 3.0 without changing significantly the accord with the experimental data. So in the calculation of the 
average polarization, we will consider the parameter f3 inside this range. One must remark that for heavier targets, 
the value of the parameters may change, as it can be seen in Fig. where /?=4.0 is reasonable for p — Ar collisions 
(mAr=39.95) and 2.0 for pp colhsions. 

The transverse rapidity parameter Pt will be assumed to be (S^j f3t = 6.6 for 800 Gev, 7.5 for 400 GeV and 8.0 for 
330 GeV. 




FIG. 12: dN/dy for 200 GeV pp and pAr collisions [io 



The reactions Ri to be considered to calculate the A polarization with eq. H23() are 

TT+A TT+A 
TT^A TT^A 

TT^A ^ ttOA . 



(30) 
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To calculate the S polarization, we consider six reactions (elastic and charge exchange), that are the charge conjugate 
of the ones shown in Fig. 
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It is important to remark that even at low energies, the exchange reactions ttA ^ ttE may occur. However, observing 
the fact that the resonance diagrams dominate the cross sections in this region, and the facts that the A* resonances 
may not be formed in ttA interactions, and that (/ettS* << gK-nT.' (considering E*(1385)) we expect that the amplitude 
of the ttA ttE reactions are not larger then 20% [23| of the ones for the ttE — > ttE, so, in a first approximation this 
kind of reactions may be neglected. 

In the case, four reactions are considered. 
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IV. RESULTS 



polarizations that are shown in 



Integrating eq. H23|) with the Monte Carlo method, we obtained the A, E and a.'^ 

Fig. cm 

Fig. 1131 shows the A as a function of pt, produced from a 400 GeV proton beam and for the production angles 
9l—7.5 and 20 mrad compared with the data from Q We can see that the model gives Pj- ~0 in the considered 
Pt region, what is consistent with the experimental results. This fact may be justified observing Fig. 01 in an average 
calculation the small contributions of the microscopic processes are almost totally canceled. 
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FIG. 13: A polarization compared with the experimental data of Q and (proton beam). 

In Fig. El the calculated polarization is now compared with the data obtained with a 330 GeV E^ beam and 
the prediction P-^ ~0 is still consistent with the data. 

The calculated S"*^ polarization can be seen in Fig. 1151 The resultant polarization is negative and with /? ^2, the 
results are in accord with the experimental data from the E761 experiment pTj |. where 9^=2.4 mrad and p^^^SOO 
GeV. 
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FIG. 14: A Polarization for 13 — 2, compared with the WA89 experimental data P|, obtained with a 330 GeV E beam and at 
< XF >=0.11. 
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FIG. 15: polarization, for some values of /3 compared with the E761 data flo| (production angle of 2.4 mrad and incoming 
beam of 800 GeV). 
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FIG. 16: S polarization compared with the data of [ll| (for < xf >=0.5). 
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FIG. 17: A polarization, for /3 = 2 and some production angles 




FIG. 18: H polarization, for /? = 2 and some production angles 
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FIG. 19: E polarization for /3 = 3 and some production angles 
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The resultant S polarization is positive, but with magnitude smaller then the one of the experimental data |ll| 
(Fig. I15|l . Considering that the experiment is for p-Cu collisions one must expect a value of /3 between 2.0 and 4.0 
(see Fig. (fT^ l. The production angle 9l is between 0.9 mrad and 2.9 mrad. and we made the calculations in this 
region, considering the average value of in the experiment, that is 0.5 |ll| . In this case, we expect that some 
effects must be included in order to improve the accord with the experimental data. 

Fig. 1 171 191 shows the A, E and polarization dependence with the production angle 9l- Observing the figures, 
we can conclude that the antihyperon polarization is a small angle effect, increasing 9^ the polarization vanishes. The 
model predicts Py ''^0 for 9^ >50 mrad. This fact explains why the antihyperon polarization was found experimentally 
in angles as small as few mrad. 

V. DISCUSSION OF THE RESULTS 

In this paper, we proposed a model in order to explain the antihyperon polarization data. It is an indirect mechanism, 
the antihyperon is produced unpolarized, and then becomes polarized by the final-state interactions. This model 
supplements the direct mechanisms, that describe the hyperon polarization, and are based on the leading particle 
effect 12J- 15] . We expect that the proposed mechanism also produces hyperons, but when compared with the direct 
mechanisms, it will be just a correction. 

The final A and S polarizations agree quite well with the experimental data. The E polarization is also in accord 
with the data, but in order to improve the agreement, some effects that we neglected during the calculation may be 
taken into account. The ttA ttS, for example were not included, and they may contribute to include the final S 
polarization. This kind of reaction will be included in future works. The pYY' coupling is another object that must 
be considered, the current value for these couplings is far from being well established yet. 

As we can see, the model shows that the Y polarization is a small angle effect, our prediction is that Py is very 
small for 9l > 50 mrad. Another result is that Py does not depends much on the incident particle {p, E ...). It may 
only change when the parametrization of da /da changes, what would occur specially in AA collisions. 
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APPENDIX A: BASIC FORMALISM 

In this paper p and p' are the initial and final hyperon 4-niomenta, k and k' are the initial and final pion 4-niomenta, 
so the Mandelstam variables are 

s = (:p + kf {p' + k'f (Ai) 

t = [p-p'f = {k-k'f (A2) 

u = {p' -kf = {p-k'f . (A3) 

With these variables, we can define 

V = (A4) 

9m'^ — f 

"« = ^ («' 

- - . (^«) 

where m, rrir and m^^ are, respectively, the hyperon mass, the resonance mass and the pion mass. The scattering 
amplitude for an isospin / state is 

Ti = + itLEle^^ u{p) , (A7) 

where Ai and Bj are calculated using the Feynman diagrams. So the scattering matrix is 



with 



m=^-^^[Aj + {^s-m)Bj] , (A9) 
fi{0)=^-^^[-Aj + {V~s + m)Bj] , (AlO) 



where E is the hyperon energy. The partial-wave decomposition is done with 

1 

ai± = -J^ [Pi{x)fi{x) + Pj±i(x)/2(x)] dx . (All) 
In our calculation (tree level) ai± is real. With the unitarization, as explained in Section III, we obtain 

These complex amplitudes are used to calculate 



/W = + + lai-] Pi{x) , (A13) 

1=0 

oo 

g{e)=iY,[ai+-ai-]Pl'\x) . (A14) 



We have, then, in the center-of-mass frame, 

da 
dQ 



=I/I' + |5|', (A15) 

f = ^ [l/P + , (A16) 

at =47r^[(/ + l)h+|2 + /|a;_|2] . (A18) 
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The YttY' coupling constants may be obtained considering the SU(3) relations |4j| 

gTvNN — 9 

2 

9T.TTT. = 2g(l - a) 

gs.s = ,9(1 -2a) , (A19) 

considering the current value (7^ArAr=13.40 and gAr^'s — 12.92, obtained in |42l | from the study of the hypernuclei, we 
have a=0.83, and the other coupling constants are determined. Similar relations hold for the Y pY' couplings. The 
YttY* couplings are calculated with the procedure described in p^ . 
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